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Abstract 

The adsorption of 4 He inside and outside a single fullerene Ceo is studied. A 

physisorption potential is proposed. The energetics and structural features of 

C6o- 4 He7v clusters are investigated. Particular attention is paid to the growth 

of the highly pronounced layered density profile. The evolution towards bulk 

liquid and surface thickness at the free interface are discussed. 
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I. INTRODUCTION 



Since the discovery of the fullerenes and nanotubes a large amount of work has been 
devoted to study properties of such systems. 1 In particular, in the case of nanotubes the 
adsorption of quantum fluids has been thoroughly investigated. 2 The physisorption poten- 
tial for such systems has been modeled several years ago by Stan and Cole. 3 Calculations 
performed for 4 He adsorbed into carbon nanotubes have shown interesting features. 4 ' 5 It was 
found that inside tubes of radius ^ 7A, besides the formation of a shell phase located at 
about 3 A from the wall, a quasi-one-dimensional structure along the cylinder axis is also 
developed. These sorts of quasi-one-dimensional structures are of theoretical interest. 6,7 On 
the other hand, the shell phase located close to the wall exhibits a surface density in agree- 
ment with the experimental value 8,9 obtained for monolayer completion density on graphite 
(see Table I). Moreover, adsorption in bundles of nanotubes have been also studied showing 
a rich pattern of phase transitions. 10 ' 11 

On the other hand, hitherto there is no study of the adsorption of 4 He onto C 60 in 
the literature. One expects that the cloud of 4 He surrounding this fullerene will exhibit 
a shell structure of the type found for helium clusters doped with atomic or molecular 
impurities. 12-18 It is of interest to determine how pronounced will be the layering as well 
as to estimate the shells' surface densities in order to make a comparison with results for 
other geometries mentioned above. By looking at the behavior of these C6o- 4 Hejv clusters 
for increasing number of helium atoms, N, one may search whether growth instabilities 
like those found in the case of planar graphite substrates 19 are also present in this case. 
Furthermore, one may search for what is going on when the clusters approach bulk liquid. 
Other issue to be examined is the liquid-vacuum interface thickness. There is a current 
interest in this observable 20 and experimental data for spherical and planar helium systems 
are available. 21,22 Moreover, since the field is open for exploration new unexpected features 
may appear. 

The paper is organized in the following way. In Sec. II we outline the model and construct 
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the physisorption potential for the interaction between Ceo and 4 He atoms. Numerical results 
are presented in Sec. Ill together with the discussion of the whole picture exhibited by the 
systems. Section IV is devoted to final remarks. 

II. THE MODEL 

A. Theoretical Framework 

The adsorption properties at T = K may be studied by analyzing the grand thermo- 
dynamic potential 

tt = E gs -pN, (2.1) 

where E gs is the ground-state energy, \i the chemical potential, and N the number of particles 
of the adsorbed fluid. In the present work we adopted a density functional (DF) approach, 
which has proven to be a successful tool for treating this kind of quantum many-body 
problems. In such a theory the ground-state energy of an interacting iV-body system of 4 He 
atoms, confined by an adsorbate-substrate potential C/ su b(r), may be written as 

+ J drp(r)U snh (r) , (2.2) 

where p(r) is the one-body density and e sc (r) the self-correlation energy per particle. The 
density profile p(r) is determined from the Euler-Lagrange (EL) equation derived from the 
condition 

5il _5{E gs [p,Vp]-p_Nl = Q _ 



5p(r) 5p(r) 

In the present work, we adopted for 4 He the non-local density functional (NLDF) pro- 
posed by the Orsay-Trento (OT) collaboration, 23 with only one change which consist in 
neglecting the non-local gradient correction to the kinetic energy term as it was done by 
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Mayol et al. 2A Although this kind of functionals is known, we present a summary for the 
sake of completeness. The self-correlation energy is expressed as 



e sc (r) =\J dv,p{v,) Vf T (l r-r/ |) + y [p(r) 



(2.4) 



with c/ 4 = -2.411857 x 10 4 KA 6 and c// 4 = 1.858496 x 10 6 KA 9 . The two-body interaction, 
VJ OT (] r — r/ |), was taken as the 4 He- 4 He Lennard- Jones (LJ) potential screened in a simple 
way at distances shorter than a characteristic distance /iot 

12 / \6" 



4e 



LJ 



if r > /i 



(2.5) 







if r < to 



OT 



with the standard de Boer and Michels parameters, 25 namely, well depth e^j = 10.22 K and 
hard core radius cr LJ = 2.556 A, in addition, the screening distance tooT = 2.190 323 A. 

Quantity p(r) is the "coarse-grained density" defined as the straight average of p(r) over 
a sphere centered at r and with a radius equal to the screening distance /iqt 



p(r) = J dr/p(r/)W(\ r - r/ |) , 



(2.6) 



where W(| r — r/ |) is the normalized step function 



W(| r-r/ |) = < 



(2.7) 



if I r - r/ I > h. 



OT 



B. The physisorption potential 

The fullerene C 60 is to a very good approximation a spherical system. 1 For constructing 
the potential due to this carbon structure we followed the procedure already applied in the 
case of a single nanotube. 3 So, we suppose that an adsorbed 4 He atom located at r, measured 
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from the Ceo center, interacts with a single substrate C atom of the spherical shell at rf via 
an isotropic LJ pair potential with standard parameters £cHe = 16.24 K and ctchc = 2.74 A 
taken from Ref. 9. Next, we consider a spherically averaged potential by assuming that the 
C atoms are uniformly distributed over the spherical surface of radius i?f u ii. Under these 
conditions, the total effect of all carbon atoms is given by 



U sub(r) = 8 IT £ C Kc 6 S Rf ull 

\ 12 / 

0"CHc \ / 0"CHc 



X 



P7Y 

Jo 



6, 



sin 9 d9. 



(2.8) 



| r — ri | I \ \ r — ri \ 

Here, 6 S = N c / (47ri?f ull ) is the uniform surface density of C atoms. After some straightfor 
ward algebra the expression for t/ su b( r ) inside the fullerene can be cast into the form 



U suh (r) = 2N c e C R- c 



x 



V R J "^"> V R Ms(Z/) 

V -Kfull / V -Kfull / 



V -Kfull / 

where M n {y) stands for the integral 



(2.9) 



sin (p dip 



(1 + v 2 - 2v cosy?) n 
1 dx 



i (1 + v 2 + 2vxY 



1 



(2.10) 



2{n - \)v L(l - u) 2 ^) (1 + z/) 2 ^ 1 ) . 

with v = r/i?f u n. On the other hand, the adsorption potential outside of the fullerene 
becomes 



U suh (r) = 2N c e CUc 



x 



M 6 (i/)-(^*) .U :! (>/) 



with v = Rhiii/r. Since in the limit v — > holds 



M n {y -> 0) = 2 , 



(2.11) 



(2.12) 



one gets the correct result at the C 60 center 



U suh (r = 0) = N c Ae LJ 



/ CTCHe V _ / tTCHe V 

V -Rfuii / V .Rfuii / 



(2.13) 



and also the appropriate asymptotic behavior very far away from the fullerene 



U suh (r » i? fu n) = N c 4 Elj 



0"CHc 



12 



/ ^CHe 



V r 



(2.14) 



We expect that the resulting potential would give a reliable description of the main features 
of the systems. 



III. RESULTS AND ANALYSIS 

In the case of a spherical symmetry the variation of Eq. (2.3) leads to the following 
Hartree like equation for the square root of the one-body helium density 

h 2 ( d 2 



n~ I cr 2 d \ r— 
~2^ [d^ + rd^l V P(? 



+ 



V H {r) + U suh (r) 



(3.1) 



which also determines \i. Here Vff(r) is a Hartree mean-field potential given by the first 
functional derivative of the total correlation energy E sc [p] 

SE sc [p] 5 



V H (r) = 



J drl p{vl)e sc {rl) . 



(3.2) 



8p(r) 5p(r) 

The expression derived for the spherically symmetric Vn(r) is given in the Appendix. Equa 
tion (3.1) was solved at a fixed number of helium atoms 



N = Att / r 2 dr p(r) 
Jo 



(3.3) 



The calculations has been carried out for i? full = 3.53 A (cf. Ref. 1) that for N c = 60 leads 
to 6, = 0.38 A - ' in agreement with the surface density used for nanotubes. As shown in 
Fig. 1 the difference between the inner and outer potential depths is dramatic. 
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A. Inner Atom 



Since the radius is very small there is enough room for only one 4 He atom inside the 
fullerene. The solution for such a system is obtained from Eq. (3.1) by using U su b(r) given 
by Eq. (2.9) and setting to zero the correlation energy. So, the problem is reduced to a 
Schrodinger equation. 

Figure 2 shows the energies and wave functions for the lowest two eigenstates together 
with the inner potential. A remarkable feature of these results is that inside of the fullerene 
a helium atom is extremely strongly bound with a ground-state energy E gs — E — —592 K. 

B. The C6o- 4 Hejv clusters 

Adsorption on the external side of the fullerene has been studied by solving the complete 
Eq. (3.1) for systems up to N = 1600 4 He atoms. In this approach we neglect the zero-point 
motion of the fullerene and treat it as a body with infinite mass. It has been established that 
this is a reasonable description in the case of helium clusters doped with sulfur hexafluoride 
molecules (SF 6 ). 16-18 One expects that this approach be even better for C 60 because it is a 
heavier and stronger attractor. 

The energy per particle, e = E gs /N, and the chemical potential, /i, are plotted in Fig. 3 
as a function of N. The size of these results can be understood by taking into account that 
the adsorption potential outside presents a well depth of about 100 K, being much smaller 
than inside as shown in Fig. 1. Examples of density profiles are displayed in Fig. 4. Small 
C 6 o- 4 Hejv clusters exhibit a typical layered structure. Two well defined layers are developed 
at N ~ 160. For larger systems begins the transition towards bulk liquid. 

Let us first focus the attention on small systems. After a small region where Q > (for 
N < 20) begins the growth of a stable monolayer film. For N > 60 starts the formation of 
the second layer and /i exhibits a kink. At this point one must be careful because a growth 
instability similar to that found for films adsorbed on strongly attractive planar substrates 19 
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might appear. This kind of instability would cause a jump of N in the adsorption isotherm. 
We have checked the possible presence of such a feature by increasing N in steps of AN = 1. 
It was found that fi is always an increasing function of iV as shown in Fig. 5. In addition, we 
analyzed whether the unstable regime is close. A transition to instabilities may be induced 
by increasing the strength of the C 60 - 4 He interaction. In turn, this may be achieved by either 
changing the potential parameters or going to fullerenes with more carbon atoms, i.e., larger 
radii. Since we would like to stick Ceo, the former way for enhancing the attraction of the 
potential was selected. Calculations have been carried out for <t C Hc = 2.98 A, this value has 
been widely used in the past for the graphite-helium interaction (see discussion in Ref. 26) 
and it increases about 10% the well depth. The new plateau corresponding to the growth of 
the second layer is shown in Fig. 5. From this drawing one concludes that in this case \i is 
also an increasing function of N indicating stability of the solutions. So, we can state that 
the solutions for <tchc = 2.74 A do not undergo growth instability and, in addition, that the 
unstable regime is not close to the adopted scenario. 

We shall now examine the surface density of the first two layers. Experimentally it has 
been determined that for a planar substrate of graphite the adsorbed superfluid 4 He presents 
two well defined layers with coverages n s i = 0.115 A~ 2 and n s2 = 0.093 A" 2 for the first and 
second layer, respectively. 8 ' 9 In the present work, surface densities were evaluated according 
to 

N- 

where < Ri > is the location of the i-layer and JVj its number of particles. Looking at Fig. 
4 one realizes that the first layer is well defined, it is centered at < R\ >= 6.41 A and 
a simple integration yields Ni = 57 leading to n sl = 0.110 A -2 . Although the position 
of the second layer may be easily obtained, < R 2 >= 9.35 A, the evaluation of N 2 is not 
straightforward because the right-hand side region of the peak is not clearly exhibited. By 
taking into account the superposition of tails of the second and third peaks we estimated 
N 2 = 90 which gives n s2 = 0.082 A -2 . In order to facilitate the comparison all the results 



for the surface density mentioned along the paper are collected in Table I. A glance at this 
table indicates a fair agreement between calculated and experimental data, although the 
present value of n s2 is slightly smaller that the experimental one, presumably due to the 
curvature of Ceo- 

For large A the systems evolve towards bulk liquid. This is merely an expected behavior. 
The real theoretical prediction is how rapidly the clusters approach that limit when A is 
increased. In addition, the comparison with other doped clusters becomes of interest. In Fig. 
6 we plotted e and /j, as a function of A -1 / 3 , which is the appropriate expansion parameter 
for spherical drops. To fix ideas we recall that e and \i for free spherical 4 He clusters obey 
the following asymptotic laws: 

e As = ^ = e B + (^j V3 ^ iV~ 1/3 , (3.5) 

and 

^=^w) A r eB+2 (S} 1,3 ^ N ~ m - (3 - 6) 

Here e B = /i = -7.15 K, p = 0.021836 A" 3 , and = 0.274 KA- 2 (see Refs. 27-31). 
These straight lines are included in Fig. 6. The results displayed in this figure indicate that 
fj, reaches much more rapidly the asymptotic trend than e. At about A ~ 500 the chemical 
potential lies, in practice, on the straight line given by Eq. (3.6). At this value of A the 
slope of /i changes indicating a phase transition to a metastable regime. The behavior of e 
is biased by the very strongly attractive potential. The dashed curve indicates how e might 
attain e B . This curve was obtained by fitting data for A > 500 to 

e = e B + (^j V3 aoo A" 1 / 3 + a c + a A" 1 , (3.7) 

which yielded a c = 0.2 K and a = —3316.4 K. The large negative value of a is caused 
by the strong attraction exerted by C 60 . Furthermore, it is important to notice that these 
coefficients also explain the behavior of /i. In this region of A the chemical potential obeys 
Eq. (3.6) due to the facts that a c is small and the huge a does not contribute because 
H = dE gs /dN. 
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The trend of e towards the asymptotic law given by Eq. (3.5) is noticeable delayed in 
comparison with diffusion Monte Carlo (DMC) results 12 for SF 6 - 4 Hejv clusters. It should 
be mentioned that these DMC values were reproduced fairly well by NLDF evaluations (see 
Fig. 1 in Ref. 18) as well as by hypernetted chain calculations (cf. Fig. 12 in Ref. 15). 
Perhaps it is worthwhile to notice that the latter figure also indicates that helium doped 
with atomic impurities like Ne, Ar, Kr, and Xe reaches even more rapidly the asymptotic 
behavior for e than the SF 6 - 4 He,/v clusters do. On the other hand, from a glance at Fig. 3 
of Ref. 18 one could conjecture that the collective energy of the dipole mode for C 60 - 4 He7v 
clusters with N even a few hundreds bigger than 1000 would still not approach zero. This 
means that the instability scenario for the dopant would not be reached for N = 1600 yet. 

Let us briefly refer to the entrance to the bulk regime. From the fit to Eq. (3.7) it 
is possible to determine that the crossing between e and p takes place at iV -1 / 3 ~ 0.041, 
see also Fig. 6. This indicates that at iV ~ 14500 the bulk could be reached. It would 
be interesting to find out whether the collective energy of the dipole mode for such large 
systems is still positive. 

Finally, we shall focus our attention on the width of the liquid-vacuum interface. At this 
interface, the density profile of 4 He changes continuously from liquid density p e to zero over 
a distance of some angstroms. The width W of a surface is defined as the distance in which 
the density decreases from 0.9pc to O.lp^. For consistency with other studies we adopted 
pe = p . The analysis was centered on systems with N > 600 because for such clusters the 
oscillations of p(r) close to the free interface are already damped enough. The present results 
for W are compared in Fig. 7 with values calculated for free 4 He droplets 20 and experimental 
data. 21 ' 22 The results obtained in the measurement of large droplets published in Ref. 21 
are represented by the reported mean thickness W = 6.4 ± 1.3 A. In the case of Ref. 22 
only the two data included by the authors in the abstract are shown. Since the latter values 
correspond to rather broad planar films, we plotted them schematically close to origin for 
the abscissa. The overall agreement exhibited in Fig. 7 is good. 
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IV. FINAL REMARKS 



The main results of this investigation about C 60 - 4 Hejv clusters may be summarized in 
the following way. Only one 4 He atom may be put inside of this fullerene and its binding 
energy is very large. The energetics of the adsorption on the external side was studied and 
the structure of the films was determined. 

The evolution of the systems for increasing N does not exhibit any growth instability 
after the completion of the first layer and, in addition, the unstable scenario is not close. 
The surface density of the first two layers are in fair agreement with experimental data for 
planar graphite. The analysis of the transition towards bulk liquid as a function of N shows 
that fj, attains more rapidly the asymptotic behavior than e. Results obtained from a fit of e 
account for this feature. The C 60 - 4 HeAr clusters approach the bulk regime much slower than 
helium clusters doped with a SF 6 molecule. 

The surface thickness at the liquid-vacuum interface obtained for rather large clusters 
matches reasonably well with theoretical evaluations for free helium spheres and experimen- 
tal data for droplets and films. 

For the future work it remains the study of excitations, which may give information on 
the location of C 60 . Furthermore, this issue is also a challenge for experimental work as 
well as Monte Carlo simulations due to the fact that such investigations may also shed some 
light on the importance of the zero-point motion of this heavy dopant. Moreover, it would 
be exciting to look for a possible levitation of Ceo in bulk helium. 
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APPENDIX: HARTREE MEAN-FIELD POTENTIAL 



In this Appendix we compile the relevant expressions derived in the adopted NLDF 
approach. For a spherical geometry the definition of the Hartree mean-field of Eq. (3.2) 
leads to 

V H (r) = J dvlp{rt) Vf T (| r - r/ |)) 
+ y[p»] 2 +^ i [p»] 3 

+ J dr/p(rf){c/ 4 p(rf) + c// 4 [p{ri) ] 2 } W(| r - r/ |) . 

(Al) 

Let us first provide the expressions of the contributions involving the "coarse-grained den- 
sity" p(r), i.e., 

p(r) = J dvtp{rt) W( | r - r/ |) , (A2) 

and 

p v (r) = Jdrfp(rf)[p(rf)] n }V(\r-rf\). (A3) 

Both these integrals may be cast into the form 

H(r) = —^ r - f dvlK{ri) 9(/i T - | r - r/ |) . (A4) 
47r/i OT J 

After introducing spherical coordinates and taking into account that the step function is 
symmetric in the azimuthal angle ip, the integration over this variable yields 



3 /■r/max fSmax 

K(r) = — — / r! 2 dr!K(r!) / sin^rf^ 

2 h QT Jrf min Jo 

x &[h 2 OT -r 2 -rf 2 + 2rrf cos6] . (A5) 



Since we shall treat helium systems located at r > /i t, the limits for the radial integral 
become r/ min = r — hor and r/ max = r + /iot- Furthermore, the step function sets 

cosfl max = r2+ I /2 ~^ T . (A6) 
Zrr! 
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Hence, one gets : 
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H(r) 



r+hoT 



rl driK(ri) 



1 - 



r — rl 



h 



OT 



(A7) 



4rh T Jr-hoT 

On the other hand, the integration of the screened LJ potential contributing to Eq. (Al) for 
r > hoT leads to 



V L H JScr (r) = J drlp(rl) Vf T (| r - rl |)) 



47re L j alj f r~ho 



r 

x 



r rr—hoT roo -\ 
I JO Jr+h OT J 



+ 
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1 

~ 2 
4 ^lj c^lj 



+ / \rldrlp(rl) 

10 / , \ 10i 
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(A8) 



It is worthwhile to notice that we expressed these quantities in a rather simple compact 
form similar to Eq. (2.16) written in Ref. 32 for planar films instead of adopting the expan- 
sion in terms of Legendre polynomials proposed in Ref. 13. The latter procedure is more 
appropriate for studying excitations of a given system. 
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TABLES 

TABLE I. Surface density of 4 He layers adsorbed on planar graphite compared with results for 
helium shells confined into a cylindrical carbon nanotube and adsorbed onto a spherical Cqo. 
Geometry Method n sl [A~ 2 ] n s2 [A~ 2 ] Ref. 

Planar Experiment 0.115 0.093 8, 9 

Cylindrical Theory 0.113 4 

Spherical Theory 0.110 0.082 PW° 

a PW stands for present work. 
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FIGURES 

FIG. 1. The left-most curve is the potential for a 4 He atom inside the fullerene evaluated with 
Eq. (2.9). The right-most curve is the potential for a 4 He atom outside of the fullerene evaluated 
with Eq. (2.11). Notice that the surface of Ceo is located at r = 3.53 A. 

FIG. 2. Energies and wave functions of the two lowest eigenstates for a single 4 He atom inside 
Ceo- Quantities ip n (r) are in arbitrary scale. The bold curve is the potential. 

FIG. 3. Energy per particle and chemical potential for C6o- 4 HeAr clusters as a function of the 
number of helium atoms. 

FIG. 4. Density profiles for a series of C 60 - 4 He Ar clusters. Results for N = 50, 100, 150, 200, 
300, 400, up to 1600 are displayed. 

FIG. 5. Chemical potential as a function of N. The plateau corresponding to the growth of 
the second layer is shown for two values of <7cHe (see text) . 

FIG. 6. Energy per particle (filled circles) and chemical potential (empty circles) as a function 
of iV -1 / 3 . Asymptotic laws for e and fi are displayed. The dashed curve indicates the trend of e 
towards the asymptotic ee = no- Stars are DMC values for the energy per particle of SF6- 4 He7v 
clusters taken from Ref. 12, the dot-dashed curve stands only to guide the eye. 

FIG. 7. Thickness of the free surface as a function of iV -1 / 3 . Triangles are data for C6o- 4 He7v 
clusters with N > 600. Full circles are results for free 4 He droplets from Ref. 20 and the solid curve 
indicates the trend of these data. The "otimes" ((g)) and the "oplus" (©) are experimental values 
from Refs. 21 and 22, respectively. 
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